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Cosmic expansion and growth histories in Gahleon scalar-tensor models of dark energy 
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We study models of late-time cosmic acceleration in terms of scalar-tensor theories generalized to 
include a certain class of non-linear derivative interaction of the scalar field. The non-linear effect 
suppress the scalar-mediated force at short distances to pass solar-system tests of gravity. It is found 
that the expansion history until today is almost indistinguishable from that of the ACDM model 
or some (phantom) dark energy models, but the fate of the universe depends clearly on the model 
parameter. The growth index of matter density perturbations is computed to show that its past 
asymptotic value is given by 9/16, while the value today is as small as 0.4. 
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I. INTRODUCTION 

Since the discovery of cosmic acceleration various pos- 
sibilities have been explored to account for this mys- 
tery. Probably the most conservative possibility is 
that the accelerated expansion arises from a cosmolog- 
ical constant or some unknown dynamical field (e.g., a 
quintessence) [1]. Instead of introducing a new compo- 
nent in "T^,y" one may alternatively consider that dark 
energy is geometric, i.e., modification to general relativity 
(GR) at long distances is responsible for cosmic accelera- 
tion. The latter class of models includes scalar-tensor 
theories, f{R) gravity [2], and the Dvali-Gabadadze- 
Porrati (DGP) braneworld [3, 4]. 

A difficulty in modifying GR at long distances lies, 
however, on short distance scales; a new scalar gravi- 
tational degree of freedom, which commonly appears in 
geometric dark energy models, must be tamed carefully 
in order to pass the stringent tests of gravity in the so- 
lar system [5, 6]. Two main ways of hiding the scalar 
degree of freedom are known to exist. The first one is 
to make the scalar effectively massive in the vicinity of 
matter. This is called the chameleon mechanism [7] and 
is utilized in viable f{R) models [8]. The second way is 
decoupling the scalar from matter in the vicinity of the 
matter sources. This is known as the Vainshtein mecha- 
nism [9], and the DGP braneworld offers a nice example 
that implements it by a non-linear self-interaction in the 
kinetic term. See also a recent paper [10] for yet another 
way of suppressing the scalar-mediated force. 

Inspired by the DGP braneworld, a class of scalar- 
tensor theories of gravity has been explored that enjoys 
self-screening of the scalar-mediated force in the vicinity 
of matter due to non-linear derivative interactions. Such 
a scalar degree of freedom is called the Galileon field in 
the original proposal [11], because the equation of motion 
for the scalar 6 is invariant under d,,<ti d,,4> + b,, on 



Minkowski spacetime. The Galileon scalar-tensor theo- 
ries have been covariantized and a general form of the La- 
grangian in curved spacetime has been considered in [12] ^ 
(see also [13, 14]). The resultant equations of motion 
have only second derivatives of the fields, and hence the 
Galileon theories are in some sense similar to Lovelock 
gravity [15]. Cosmology based on a Galileon field has 
been studied in Refs. [16-18]. 

In this paper, we study aspects of Galileon cosmology 
in terms of a certain wider class of the Lagrangian than 
previously studied. The action we consider is of the form 



S 



(1) 



where (j) is the Galileon field and the coupling A(0) has 
dimension of length. Matter (represented by the La- 
grangian Cm) is universally coupled to gravity in the 
Jordan frame. The gravitational field equations derived 
from (1) are 
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and the equation of motion for the Galileon field is given 
by 



R + uj 
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^ Upon covariantization the original symmetry under d^ip 
bfj, is lost. However, covariant Galileon Lagrangians are uniquely 
determined by requiring that the equations of motion have only 
second derivatives of the fields. Thus, the name "Galileon" may 
be misleading and inappropriate. Nevertheless, we use the name 
in this paper because no other name is as good as "Galileon." 
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2A 



= 0, (3) 



where a prime denotes differentiation with respect to (p. 
The above modified gravity theory may be regarded as 
the generalization of Brans-Dicke gravity, but the cubic 
derivative interaction, [A^/0^]n0(V(/')^, plays a key role 
in manifesting the Vainshtein mechanism. The Vain- 
shtein radius, below which the Galileon-mediated force 
is screened, is evaluated as ry ~ [rgX'^/{l + 2a;/3)^]^/^, 
where rg is the Schwarzschild radius of the matter source. 

In what follows we shall consider the particular case 
with u! = const and A^ oc where a may be positive 
and negative. The self-accelerating de Sitter universe 
of [17] corresponds to the subclass a = 0. We are go- 
ing to study the modified gravity theory with general a 
as a possible alternative to a cosmological constant or 
dynamical dark energy. Accelerating cosmologies from 
a scalar field non-minimally coupled to gravity without 
the non-linear derivative interaction have been studied 
extensively in [19, 20]. 

Since the background expansion history in successful 
modified gravity is by definition almost identical to that 
of the standard ACDM model or other dynamical dark 
energy models, it is important to study the growth his- 
tory of perturbations as a tool to distinguish modified 
gravity from models with a cosmological constant or dark 
energy. We compute the growth index of matter density 
perturbations [21, 22] in Galileon cosmology, which is 
known to be a powerful discriminant among models of 
cosmic acceleration [23-25]. 

The paper is organized as follows. In the next section 
we study the background evolution of Galileon cosmol- 
ogy derived from the above action. Then, in Sec. Ill, we 
discuss the growth density perturbations in Galileon cos- 
mology, paying particular attention to the growth index. 
The final section is devoted to conclusions. 



II. COSMIC EXPANSION IN GALILEON 
THEORIES 



Let us start with investigating the background evo- 
lution of Galileon cosmology in the presence of cold 
dark matter. It is convenient to define a new field by 
(j) = (Mpj/2) cxp((t). For a flat Friedmann-Robertson- 
Walker universe, ds^ = — dt^-|-a^(f)dx^, the gravitational 



field equations are 
3H 



2 _ P^e--_sH&+'^&^ + 3X^Ha^ 
2 



to 



:-2" (A^),cr.. 



(4) 



(5) 



and the Galileon field equation is 

6 (2H'^ + h'^ - {2a + a'^ + 6H&) - uj^^a'^ 



6A2 3H^& + H& + 2Ha + ) & 



+4 [(A^),, - 2A2] a'a + [(A^),,, - (A^),, - 2X^] a' 
= 0. (6) 

Here, an ovcrdot represents differentiation with respect 
to t and H{t) :— d/a is the Hubble expansion rate. The 
energy density of cold dark matter, pm, is conserved in 
the Jordan frame, pm + S-ffpm = 0. The Friedmann equa- 
tion (4) can be written in the form 

3J^' = ^(Pm + Peff), (7) 

where the effective dark energy density Peff is defined as 



Peff 



-3H& + + 3X^H&^ 



{X% - 2X-) a' 



3ij2(i_e'^). (8) 



We assume that ^,A^„,A^^g. < O(A^), where is 
the Hubble rate at present. At early times, H~'^ <C A^, 
Eq. (6) is approximately solved to give 



X'u 



2H^ +H 
3H^ + H 



1 

3' 



(9) 



where we used H « 2/3t. (To be precise, this is the 
consequence of Pm ^ C'(Peff)) which will be verified 
shortly.) From this we see that if the variation of A 
is small over the Hubble time scale, a is constant and 
(T w at w 2&/3H. Note that with the rescaling of Mpi 
we may set a{t = 0) = 0. Two branches of solutions 
are present, d > and & < 0, but it turns out that 
the negative a branch is plagued with ghost instabili- 
ties and only the positive & branch is healthy [16] (see 
the Appendix). We therefore focus on the case with 
(7 > in this paper. In this early-time regime, we have 
/9eff/M|i « -2Hae'' + 3H'^{1 - e") « ~AHa, so that 
the usual matter-dominated universe is indeed consis- 
tently reproduced, 3M|[7J^ ~ > 0{pcfi)- This is 
the cosmological version of the Vainshtein effect [16]. It 
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is worth noting that for ct > the effective dark en- 
ergy density is negative in the initial stage. The behav- 
ior peff ~ —AHa ^ —H/X reminds us of the modified 
Friedmann equation in the normal branch of the DGP 
model [4]. 

We then consider the late-time evolution of Galileon 
cosmology. To be more specific, we study the case with 
uj = const and 



A2 



/ 2<j) 



\ hp 

\JUpj 



(10) 



The model with cj = and a = —I was studied in 
Ref. [16]. It was shown in Ref. [17] that the model with 
u! < —4/3 and a = allows for a late-time de Sitter 
solution, mimicking the ACDM model. We generalize 
the result of [17] and show that the more general class 
of models described by (10) can give rise to the cosmic 
expansion history similar to the ACDM model. 

For general a one finds that Eq. (5) admits the follow- 
ing solution: 



JL 1 



a 
H 



(11) 



where 



* ■ uj + 2-a ' ^ 

In order for s to be real, we require u! < — (2 — Q!)(4+a)/6. 
This includes the parameter range in which the scalar 
field would be a ghost in the usual Brans-Dicke theory, 
i.e., UJ < —3/2, but the situation is drastically changed by 
the introduction of the non-linear derivative interaction 
and the model does not suffer from the pathology (see 
the Appendix). For negatively large w, we have s ^ 
0(|w|-i/2) ^ I and hence \H/H'^\ < 1, leading to a 
quasi-de-Sitter expansion. The Hubble rate is explicitly 
given hy H = {ast + C) , where C is an integration 
constant. Equation (4) with Pm = gives X^{(t)H^ — 
(2 + s)/s'^. The scale rc must be tuned to satisfy this 
relation. 

The solution (11) is our candidate for the future 
asymptotic solution. In order to work out the background 
evolution from the beginning of the matter-dominated 
stage and to see the future asymptotic behavior of the 
cosmological solution, Eq. (5) and (6) are solved numeri- 
cally with the initial conditions given by the above early- 
time solution. In Fig. 1 we show the dimensionless phys- 
ical distance, r{z) := Hq dz' / H {z') , for oj = —500 and 
various a. Here, z is the redshift and the present time 
to is determined from Pm/3A/pi-ff^|t=to = f^mOj where we 
use r^mo = 0.3 throughout the paper. We see that the 
cosmic expansion histories in the Galileon models mimic 
that of the ACDM model, and the difference is less than 
5 percent for the plotted examples. A closer look shows 
that Galileon cosmology is more similar to phantom dark 
energy models at the background level, as is obvious also 
from the behavior of the effective equation of state Wcs 
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FIG. 1: Dimensionless physical distance for a; 
various a. 
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FIG. 2: Effective equation of state u)eff as a function of z for 
Lo — —500 and various a. 



that will be investigated shortly. One can confirm that 
the initial deviation of a from the value given by Eq. (9) 
does not change the result. 

The effective equation of state, Wcs := —pcs/3Hpcs — 
1, is shown in Fig. 2 for u = —500 and various a. We 
plot the range t ^ to so that we can see the asymptotic 
behavior 1 + Wctf — > ~2H/3H^ Ki as/3. A phantom- 
like behavior is found at low z, but 1 -I- Wcs gets posi- 
tive eventually for a > 0, crossing the phantom divide. 
Therefore, the phantom-like behavior is temporary and 
the big rip singularity can be avoided in this case. In con- 
trast, 1 -|- WcB remains negative in the future for a < 0, 
leading to the singular fate of the universe, H oo. 

Since pcs oc iJ ~ 2/3t at early times, we have WcS — >■ 
— l/2as2:— ^oo. The effective dark energy density is neg- 
ative initially and evolves into pcff > at lower z, which 
implies that peff crosses zero at some z = and we have 
WcS — >■ ±(X) as 2 — >■ z^. However, this divergence is an 
artifact of the definitions of poff and WcS and no anomaly 
is seen in the evolution of H(z). The whole behavior of 
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Wcff is most similar to what is found in the normal branch 
of the DGP braneworld with the additional cosmological 
constant or a quintessence field on the brane [26-28]. 



III. DENSITY PERTURBATIONS 

A. Some preliminEiries 

In various dark energy models and modified gravity 
models the governing equations for the homogeneous 
background and the cold dark matter density perturba- 
tion 6 can often be written in the form 



6 + 2HS = 



■ Pm + Peflf, 

^6. 



2M2i 



(13) 
(14) 



Here, ^ = (_{t) represents the ratio of the effective grav- 
itational coupling "Gcff" to G[:— l/(87rMpj)], which in 
general differs from unity. ^ We shall show in the next 
subsection that the perturbation equation indeed takes 
the form of (14) in Galileon theories. Before going to the 
analysis of the specific cases of Galileon modified gravity, 
we explain how we can characterize the growth of density 
perturbations that follows from Eqs. (13) and (14). 

We are going to study the growth index 7 of density 
fluctuations. The growth index is defined via the growth 
rate / = dlnJ/dlna as 



In/ 



Inrim(ci) ' 



(15) 



where rim (a) = /Om/(3M|j_ff^). This quantity was pio- 
neered by [21, 22], and it turned out that 7 is a powerful 
and useful discriminant among different models of cos- 
mic acceleration as 7 has a different but almost constant 
value depending on the model [23-25]. For instance, 
7 ~ 0.55 for the ACDM model and 7 ~ 0.68 for the 
DGP model. The purpose of this subsection is to red- 
erive the useful formula for the asymptotic value of 7, 
7(z 00) =: 7oo [24]. 

In terms of the the growth factor, g{a) := 6/ a, Eq. (14) 
can be written as 



da2 



1 



5 + 



dlnH 
din a 



3 + 



dlnJ? 
din a 



d5 
da 

[^(a)f^m(a) 



5 = 0. (16) 



^ The cfTcctivc gravitational coupling may depend on the scale k, 
but we do not care about this point as we are not interested in 
the power spectrum of S in this paper. 



We now introduce the new variable defined by^ 

1 



y 



^m(a) 



1. 



(17) 



In another way, it is written as y = Peff/Pm- In terms of 
y Eq. (16) can be written as 



f_9 
dy 

3 + C 3 ^{y)-l 



+ 



. 5/35, 



dg 
dy 



2C2 2(2 



y 



0, 



(18) 



where ( := dlny/dlno and ( :— [(1 + t/)/(^]dln(/dlna. 

In order to focus on analytically tractable cases, we 
consider the regime in which the following two assump- 
tions are valid. The first assumption is that ( = const 
and hence C = 0. If dark energy is a fluid whose equation 
of state is w = const, then C — —3w and the assump- 
tion is valid over the whole history of the universe after 
matter domination.* If instead the accelerated expansion 
is cause by modification of gravity, ( is not constant in 
general, but at early times {y <^ 1) it is often good to 
approximate ( = const. 

The second assumption is that ^ can be approximated 
by ^ — 1 = Ay + .... We do not consider the case where 
the leading term is given by y" with n < 1, because if so 
then 6 cannot have the desired early-time behavior 5 ~ a 
(i.e., g w 1). If the leading term is with n > 1 or 
^ = 1 exactly, we neglect the effect of this term and set 
^ = 0. An interesting example for which ^ ^ is the 
DGP braneworld model [24, 29, 30]. 

Within the range of the validity of the above two as- 
sumptions, the solution to Eq. (18) is given by 



3 + C- 



1 



C(5 + 2C) 



1 + 5/(20; 

3^ 

y + -, 



-V) 



(19) 



where 




The expression in the second line of Eq. (19) is sufficient 
for our purpose. However, the expression in the first line 
might be useful if one considers the model with ^ = 1 
(i.e., A = Q) and C = const over the whole history of the 
universe after matter domination because the expression 



If the scale factor is a multiple-valued function of y, it is not 
appropriate to use y as a global time variable. This is indeed the 

case for the present Galileon model. Even so, the scale factor is 
a one-valued function of y for sufficiently small a. 
In the w = const fluid model of dark energy, the evolution of (5 
will be governed by Eq. (14) provided that the effect of clustering 
of the dark fluid is negligible. 



5 



in terms of the hypergeometric function then gives the 
exact solution [31-34]. 

Now the asymptotic value of the growth index can com- 
puted as 7cx3 = -ln(l + Cdln3/dlny)/ln(l + y)\y^o = 
—(dg/dy\y^o- Thus, we obtain [24] 



'Too 



3 + C - 3A 
5 + 2C 



B. Growth of matter density perturbations 

We now turn back to Galilcon modified gravity and 
study the growth index of density perturbations. The 
metric perturbations in the Newtonian gauge is written 
as 



(1 + 2<^)df + a^{l- 2*)dx2 



(22) 



The perturbed energy-momentum tensor of cold dark 
matter is given by 



-Spn 



ST,' 



diSq„ 



ST,' = 0. 



(23) 



and the perturbation of the Galileon field is described 
by (7 ^ a{t) + Sa{t,x). From the traceless part of the 
gravitational field equations we immediately obtain 



^ - ^ = 5(7. 



(24) 



This is in contrast to the relation one finds in general 
relativity, ^ = ^. 

We are going to investigate the subhorizon evolution 
of the comoving density perturbation, S := {5pm — 
SHSqm)/ Pm, by employing the quasi-static approxima- 
tion. The quasi-static approximation was first introduced 
in [35] in the context of quintessence cosmological mod- 
els and then was applied to various theories of modi- 
fied gravity [20, 30, 36]. The approximation amounts to 
V^X/a^ > H'^X > HX,X, where X = $,*,^o-, and 
y2 S^iQ^Q-, Although the full cosmological perturba- 
tion equations are involved [18], the quasi-static approx- 
imation simplifies the field equations to 



a2 



0, 



(25) 
(26) 



where 



So ■■ = 



2* - $ -I- uSa 

+X'^[&^<i> + 2{a + &'^ + 2Ha)Sa] , (27) 
2^~{1- A^a^) Sa. (28) 



The first equation is derived from the Galileon equation 
of motion and the second from the gravitational field 
equations. Note that the above equations are indepen- 
dent of the detailed structure of the coupling function 



A^(cr) in the sense that they do not include the terms 
(A^),^ and {X'^),aa- Using Eqs. (24)-(26) we obtain the 
modified Poisson equation. 



Pn 



2M2j 



1 



T 



S, 



(21) where 



T := 3 + 2oj + X'^ [4o- -|- 2(t^ -|- 8Ha - A^(t^ 



(29) 



(30) 



The perturbed energy-momentum conservation equa- 
tions, S {VuT^) = 0, yield 



5 + 2H5= ^<^. 



(31) 



Equations (29) and (31) are combined to give the evo- 
lution equation for the density perturbation 5 in the 
Galileon models [16, 17], which is of the form (14) with 
the identification ^ = €-"[1 + {I + X^a'^f/F]. We have 
^—1 ^ 1 initially, and then ^ increases with time, leading 
eventually to ^ — 1 = 0(1) today. This effect enhances 
the growth rate of density perturbations [17, 18]. 

As we already have the formula (21) it is easy to com- 
pute the asymptotic growth index analytically. At 
early times we have X^cr^ « 1/3, e"^ ~ 1 -I- 2a /iH, and 

« 8H/3&. Noting that y = Peff/Pm ~ —Aa/SH, we 
find ^ w 1 + 0(2/^) and dlnt//dlna w 3/2. This corre- 
sponds to C = 3/2 and A = in the formula (21). We 
thus arrive at 



7cx, = = 0.5625. 
lb 



(32) 



This is certainly different from the ACDM value, 
7oo,ACDM = 6/11 = 0.5454..., but unfortunately the dif- 
ference is only 3 percent. The value is also different from 
the DGP value, 7oo.dgp — 11/16, because of the differ- 
ent rate of time-variation in the effective gravitational 
coupling. 

Numerical solutions to Eq. (31) with Eq. (29) are also 
easy to obtain. Solving the background equations first as 
was done in the previous section and then integrating the 
perturbation equation under the initial condition 5 « a 
at t « 0, we have computed the growth index numeri- 
cally. The results are presented in Figs. 3 and 4. One 
can confirm that the asymptotic value is indeed given by 
7oo = 9/16. Since f2m(a) = 1 when peff = 0, we have 
7 — > ±00 as z — > 2;^. Similarly to the behavior of Wp[[ 
at z = z*, this is an artifact of the definition and the 
evolution of 6 is not ill-behaved. The detailed behavior 
at lower z depends upon the parameters lo and a, but 
typically, 7 is as small as ~ 0.4 due to the enhanced 
growth of density perturbations. Such a small value of 
7 due to the enhanced growth of S is also reported in 
some f{R) gravity models [37].^ In any case, the growth 



The growth index in f{R) gravity depends on the scale k, in 
contrast to the present Gahleon models. 
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FIG. 3: 7 as a function of z for a = 4 and a — —1. The 
parameters are given by a; = —500 and fimo = 0.3. 
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FIG. 4: 7 as a function of z for a — 1 and various uj. The 
parameters are given by Slmo = 0.3. 

index shows a distinct feature compared to the ACDM 
model and other modified gravity models such as self- 
accelerating DGP braneworld. 



FIG. 5: The sound speed squared (defined in the Appendix) 
as a function of z. The sound speed squared is positive, and 
can be temporarily greater than 1. 

dark energy models. We have clarified that the late-time 
asymptotics depends on a: a future big rip singular- 
ity occurs for a < 0, while the phantom-like behavior 
is only temporal and we eventually have WcS > — 1 for 
Q > 0. The result here is the generalization of the self- 
accelerating de Sitter universe found for a = in [17]. 

We have also investigated the growth history of density 
perturbations in Galileon cosmology. The asymptotic 
growth index was obtained analytically as = 9/16, 
which was confirmed numerically. The growth index to- 
day was found to be 7 ^ 0.4. The behavior of the 
growth index is thus a clear discriminant among Galileon 
scalar-tensor gravity, the ACDM model, and the DGP 
braneworld. 
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IV. CONCLUSIONS 



Appendix A: Stability analysis 



In this paper, we have studied a class of Galileon 
cosmology, i.e., cosmology in the Brans-Dicke theory 
with the non-linear derivative interaction of the form 
[\'^{<j))/(t>^]{y<f))^a<f) with A2(0) cx 0". Thanks to this 
term in the Lagrangian, the Vainshtein mechanism works 
and the scalar-mediated force is screened near matter 
sources. The Vainshtein mechanism operates in the cos- 
mological context as well, leaving the cosmic expan- 
sion history indistinguishable from the usual matter- 
dominated universe. Moreover, we have illustrated that 
modification to gravity due to the Galileon field can be 
responsible for the late-time acceleration of the universe, 
closely mimicking the ACDM model or rather phantom 



In the main text we advocate the ghost-free branch 
a > 0. The absence of ghost instabilities can be proved 
by expanding the action at second order in the fluctuation 
of the Galileon field, 6a. Noting that the cubic term 
in the action is dominant at early times, we obtain the 
following quadratic action [16]: 



(Al) 



From this we see that tr > is required in order to avoid 
ghost instabilities. (The convention here is different from 
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the one used in [16]. Our a is related to tt in [16] as 
cr -27r/Mpi.) 

Following Ref. [17] we neglect the friction term (a 5a) 

and the effect of metric perturbations on small scales to 
obtain the evolution equation for 6a: 

g{t)5a - T{t)-^5a = 0, (A2) 

where T is defined in Eq. (30) and 

g := 3 + 2uj + {UHa + 2a^ + SA^ct^) - 4{X'^)^„a^. 

(A3) 

(Full perturbation equations are found in the Appendix 
of Ref. [18].) It is easy to check that for iJ^ ^ 1/^2 
Eq. (A2) reduces to the equation of motion derived from 



the action (Al). The condition for avoiding ghost insta- 
bilities is given hy Q > 0. We also require := J^/G > 
so that the Galileon field is pertin-batively stable. Both 
are satisfied in our numerical examples presented in the 
main text thanks to the non-linear term in the La- 
grangian (see Fig 5). Note that the sound speed Cs can 
be greater than the speed of light, but that is a tem- 
porary behavior. Super luminal propagation might in- 
dicate pathology and is claimed to be problematic in 
Refs. [38-41]. However, some authors argue that propa- 
gation faster than light does not necessarily lead to prob- 
lems [42-44]. Thus, whether a theory with superluminal 
propagation is viable or not is still a controversial issue. 
We see that — >• as z -> — 1 and no particular behavior 
is found at the phantom divide crossing. 



[1] For comprehensive review, see, e.g., E. J. Copcland, 

M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15, 

1753 (2006) [arXiv:hep-th/0603057]. 
[2] For comprehensive review, sec, e.g., T. P. Sotiriou and 

V. Faraoni, arXiv: 0805. 1726 [gr-qc]; A. De Felice and 

S. Tsujikawa, arXiv: 1002.4928 [gr-qc] . 
[3] G. R. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. 

B 485, 208 (2000) [arXiv:hep-th/0005016]. 
[4] C. DeflFayet, G. R. Dvali and G. Gabadadze, Phys. Rev. 

D 65, 044023 (2002) [arXiv:astro-ph/0105068]. 
[5] C. M. Will, Living Rev. Rel. 9, 3 (2005) [arXiv: gr-qc/ 

0510072]. 

[6] J. P. Uzan, arXiv: 0908. 2243 [astro-ph.CO]. 

[7] D. F. Mota and J. D. Barrow, Phys. Lett. B 581, 
141 (2004) [arXiv:astro-ph/0306047]; J. Khoury and 
A. Weltman, Phys. Rev. Lett. 93, 171104 (2004) [arXiv: 
astro-ph/0309300]; J. Khoury and A. Weltman, Phys. 
Rev. D 69, 044026 (2004) [arXiv:astro-ph/030941l]. 

[8] W. Hu and I. Sawicki, Phys. Rev. D 76, 064004 
(2007) [arXiv: 0705. 1158 [astro-ph]]; A. A. Starobinsky, 
JETP Lett. 86, 157 (2007) [arXiv : 0706 . 2041 [astro- 
ph]]; S. A. Appleby and R. A. Battye, Phys. Lett. B 
654, 7 (2007) [arXiv: 0705. 3199 [astro-ph]]; S. Nojiri 
and S. D. Odintsov, Phys. Rev. D 77, 026007 (2008) 
[arXiv: 0710. 1738 [hep-th]]. 

[9] A. I. Vainshtein, Phys. Lett. B 39, 393 (1972). 
[10] K. Hinterbichler and J. Khoury, arXiv: 1001 .4525 [hep- 
th]. 

[11] A. Nicohs, R. Rattazzi and E. Trincherini, Phys. Rev. D 
79, 064036 (2009) [arXiv: 0811 .2197 [hep-th]]. 

[12] C. Deffayet, G. Esposito-Farese and A. Vikman, Phys. 
Rev. D 79, 084003 (2009) [arXiv: 0901 . 1314 [hcp-th]]; 
C. DefFayct, S. Deser and G. Esposito-Farese, Phys. Rev. 
D 80, 064015 (2009) [arXiv: 0906. 1967 [gr-qc]]. 

[13] C. Charmousis, S. C. Davis and J. F. Dufaux, JHEP 
0312, 029 (2003) [arXiv:hep-th/0309083]; C. Char- 
mousis, Lcct. Notes Phys. 769, 299 (2009) [arXiv:0805. 
0568 [gr-qc]]. 

[14] E. Elizalde, A. G. Zhckscnacv, S. D. Odintsov and 
I. L. Shapiro, Class. Quant. Grav. 12, 1385 (1995) 
[arXi V : hep- th/941 206 1] . 



[15] D. Lovelock, J. Math. Phys. 12, 498 (1971). 

[16] N. Chow and J. Khoury Phys. Rev. D 80, 024037 (2009) 
[arXiv: 0905. 1325 [hep-th]]. 

[17] F. P. Silva and K. Koyama, Phys. Rev. D 80, 121301 
(2009) [arXiv: 0909. 4538 [astro-ph.CO]]. 

[18] T. Kobayashi, H. Tashiro and D. Suzuki, Phys. Rev. D 
81, 063513 (2010) [arXiv: 0912. 4641 [astro-ph.CO]]. 

[19] Sec, e.g., J. P. Uzan, Phys. Rev. D 59, 123510 (1999) 
[arXiv:gr-qc/9903004]; L. Amendola, Phys. Rev. D 
60, 043501 (1999) [arXiv:astro-ph/9904120]; T. Cliiba, 
Phys. Rev. D 60, 083508 (1999) [arXiv:gr-qc/9903094]; 
F. Perrotta, C. Baccigalupi and S. Matarrese, Phys. 
Rev. D 61, 023507 (2000) [arXiv:astro-ph/9906066]; 
O. Bcrtolami and P. ,1. Martins, Phys. Rev. D 61, 
064007 (2000) [arXiv:gr-qc/9910056]; G. Esposito- 
Farese and D. Polarski, Phys. Rev. D 63, 063504 (2001) 
[arXiv:gr-qc/0009034]; E. Gunzig, A. Saa, L. Brenig, 
V. Faraoni, T. M. Rocha Filho and A. Figueiredo, 
Phys. Rev. D 63, 067301 (2001) [arXiv:gr-qc/0012085]; 
T. Chiba, Phys. Rev. D 64, 103503 (2001) [arXiv: 
astro-ph/0106550]; D. F. Torres, Phys. Rev. D 66, 
043522 (2002) [arXiv : astro-ph/0204504]; R. Gannouji, 
D. Polarski, A. Ranquet and A. A. Starobinsky, JCAP 
0609, 016 (2006) [arXiv:astro-pli/0606287]; S. Carloni, 
S. Capozziello, J. A. Leach and P. K. S. Dunsby, Class. 
Quant. Grav. 25, 035008 (2008) [arXiv:gr-qc/0701009]; 
S. Tsujikawa, K. Uddin, S. Mizuno, R. Tavakol and 
J. Yokoyama, Phys. Rev. D 77, 103009 (2008) [arXiv: 
0803.1106 [astro-ph]]; T. Chiba, M. Siino and M. Yam- 
aguchi, arXiv : 1002 . 2986 [astro-ph.CO]. 

[20] B. Boisseau, G. Esposito-Farese, D. Polarski and 
A. A. Starobinsky, Phys. Rev. Lett. 85, 2236 (2000) 
[arXiv : gr-qc/0001066] . 

[21] P. J. E. Peebles, The Large-Scale Structure of the Uni- 
verse (Princeton University Press, Princeton, New Jer- 
sey, 1980); P. J. E. Peebles, Astrophys. J. 284, 439 
(1984). 

[22] L. M. Wang and P. J. Steinhardt, Astrophys. J. 508, 483 

(1998) [arXiv:astro-ph/9804015]. 
[23] E. V. Linder, Phys. Rev. D 72, 043529 (2005) [arXiv: 

astro-ph/0507263] . 



8 



[24] E. V. Lindcr and R. N. Cahn, Astropart. Phys. 28, 481 

(2007) [arXiv:astro-ph/0701317]. 

[25] D. Polarski and R. Gannouji, Phys. Lett. B 660, 439 

(2008) [arXiv: 0710. 1510 [astro-ph]]. 

[26] V. Salini and Y. Shtanov, JCAP 0311, 014 (2003) 

[arXiv : astro-pli/0202346] . 
[27] A. Lue and G. D. Starkman, Phys. Rev. D 70, 101501 

(2004) [arXiv:astro-ph/0408246]. 
[28] L. P. Chimcnto, R. Lazkoz, R. Maartens and I. Quiros, 

JCAP 0609, 004 (2006) [arXiv:astro-ph/0605450]. 
[29] A. Lue, R. Scoccimarro and G. D. Starkman, Phys. Rev. 

D 69, 124015 (2004) [arXiv:astro-pli/0401515]. 
[30] K. Koyama and R. Maartens, JCAP 0601, 016 (2006) 

[arXiv : astro-ph/051 1634] . 
[31] D. J. Heath, Mon. Not. R. Astron. Soc. 179, 351 (1977). 
[32] S. Bildhaucr, T. Buchcrt and M. Kasai, Astron. Astro- 

phys. 263, 23 (1992). 
[33] V. Sahni and P. Coles, Phys. Rept. 262, 1 (1995) [arXiv: 

astro-pli/9505005] . 
[34] S. Lee and K. W. Ng, arXiv: 0906. 1643 [astro-ph. CO]. 
[35] A. A. Starobinsky, JETP Lett. 68, 757 (1998) [Pisma 

Zh. Eksp. Teor. Fiz. 68, 721 (1998)] [arXiv: astro-ph/ 

9810431]. 

[36] S. Tsujikawa, Phys. Rev. D 76, 023514 (2007) [arXiv: 
0705.1032 [astro-ph]]; S. Tsujikawa, K. Uddin and 



R. Tavakol, Phys. Rev. D 77, 043007 (2008) [arXiv: 

0712.0082 [astro-ph]]. 
[37] R. Gannouji. B. Moracs and D. Polarski, JCAP 0902, 

034 (2009) [arXiv: 0809. 3374 [astro-ph]]; S. Tsujikawa, 

R. Gannouji, B. Moraes and D. Polarski, Phys. Rev. 

D 80, 084044 (2009) [arXiv : 0908 . 2669 [astro-ph.CO]]; 

T. Narikawa and K. Yamamoto, arXiv : 0912 . 1445 

[astro-ph.CO]; H. Motohashi, A. A. Starobinsky and 

J. Yokoyama, arXiv: 1002. 1141 [astro-ph.CO]. 
[38] A. Adams, N. Arkani-Hamed, S. Dubovsky, A. Nico- 

Us and R. Rattazzi, JHEP 0610, 014 (2006) [arXiv: 

hep-th/0602178]. 
[39] C. Bonvin, C. Caprini and R. Durrer, Phys. Rev. Lett. 

97, 081303 (2006) [arXiv:astro-ph/0606584]. 
[40] G. Elhs, R. Maartens and M. A. H. MacCallum, Gen. 

Rel. Grav. 39, 1651 (2007) [arXiv:gr-qc/070312l]. 
[41] C. Bonvin, C. Caprini and R. Durrer, arXiv : 0706 . 1538 

[astro-ph] . 

[42] E. Babichev, V. Mukhanov and A. Vikman, JHEP 0802, 

101 (2008) [arXiv: 0708. 0561 [hcp-th]]. 
[43] V. Gorini, A. Y. Kamenshchik, U. Moschella, O. F. Pi- 

attclla and A. A. Starobinsky, JCAP 0802, 016 (2008) 

[arXiv: 071 1.4242 [astro-ph]]. 
[44] R. Geroch, arXiv: 1005. 1614 [gr-qc]. 



